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Measurements by the Belle Collaboration of the exclusive J/ip + r\ c and inclusive J/ip + 
cc+ X productions in e + e - annihilation differ substantially from theoretical predictions 
based on the nonrelativistic QCD factorization approach. In order to test if such a 
discrepancy is originated from the large perturbative corrections to the hard-scattering 
amplitude, we study inclusive production of four charm hadrons in e + e~ annihilation at 
B factories. 
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| The nonrelativistic QCD (NRQCD) factorization formalisrrP has enjoyed con- 

siderable success in describing production and decay rates of heavy quarkonia. The 
I | approach provides an infrared-safe prediction for the P-wave quarkonium decajGl 

and explains the large empirical cross sections of prompt charmonia at the Fermi- 
f~| , lab TevatronPI However, there are a few serious challenges to NRQCD. The first 

issue is the discrepancy between the NRQCD prediction and the CDF data for 
r ■ the polarization of prompt J/ip in large-pT bins!^ The measurement also confronts 

l^j , with a rec ent lattice NRQCD calculation, which supports the dominance of trans- 

' verse 

j/^EH Th e second issue is for the J/tp production in e + e annihilation at 
S-factories. The cross sections for exclusiveproduction of J/ip + ilc m e + e~ annihi- 
lation measured by the Belle CollaboratiorPB and by the Babar CollaboratiorP are 
greater than NRQCD predictions of leading-order by an order of magni- 

tude. A few proposals^ introduced to solve the problem were disfavored by a recent 
experimental analysis0The next-to-leading-order QCD corrections to the exclusive 
process have been found to be as large as 809^121 f the lea ding -order prediction!^ 
It is argued that nonperturbative corrections may be large. A * However, these cor- 
rections are not yet sufficient to explain the data. It would be interesting to see 
if such a large per turbative correction is also true in inclusive production of four 
charm hadronsp^ having a similar parton-level process. 

Another large discrepancy between the NRQCD prediction and the B- 
factory data is for the inclusive prompt J/ip + cc production cross section 
cr(e + e~ — > J/ip + cc + X). Mea sured cross sectiorP is significantly larger than 
the NRQCD predictions! 6 LLZI i n order to find if an alternative way could re- 
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Fig. 1. (a) Total cross section <r(e+e — — > cccc) at ^/s = 10.6 GeV in pb as a function of m c , 
where a = 1/137 and ct s = 0.2. (b) Differential cross section dcr/dm cc in fb/GeV with respect 
to the invariant mass m cc = mi3 of cc for e + e~ annihilation into cccc, where m c = 1.5 GeV, 
a = 1/137, and a s = 0.2. Physical range of the m cc is from 2m c to y/s — 2m c . The area under the 
curves are the integrated cross sections 97 fb. From Ref. 14. 

solve the prob lem, a proc es was recently studied within the color-evaporation 
mod el (GEM).™™ 

The CEM prediction for the cross section cr(e + e — > 
J/^ + cc + X) was reportedly smaller than the empirical value by about two 
orders of magnitude!^ A more comprehensive discussion of the CEM can be 
found in Refs l23l 1241 1251 The Belle Collaboration also measured the fraction 
R[J/ip + cc] of the process within inc lusive J/ip production. The measured value, 
R[J/ip + cc] Belle = 0.82 ± 0.15 ± 0.14j^l is significantly larger than the NRQCD 
predictions of about 15%. In the CEM, there is only a single universal nonper- 
turbative factor, which perfectly cancels in the ratio R[J/4> + cc]. Thus, the CEM 
prediction for the ratio is expected to be more reliable than the absolute values for 
the cross sections in the numerator and the denominator. Here we briefly review 
the production of four charm hadrons in e + e~ annihilation at B factories. 



Since the inclusive four charm hadron production involves the same Feynman 
diagrams for exclusive J/ip + rjc production, measuring the cross section will provide 
an important information in the estimating the size of the short-distance coefficient 
for J ftp + rj c cross section. We present er(e + e~ — * cccc + X) prediction in leading 
order of strong coupling constant, a 2 a 2 s . 

Our predictions for the inclusive four charm hadron cross sections in e + e~ an- 
nihilation at a/s =10.6 GeV depending on the charm quark mass m c is shown in 
Fig. ID^a). The cross section for e + e~ — * cccc is very sensitive to the value of m c . 
For a — 1/137, a s — 0.2, m c =1.5 GeV <r(e + e~ — * cccc) —97 fb. The cross section 



2, 2008 4:4 WSPC/INSTRUCTION FILE qcd05beijing 



Inclusive Production of Four Charm Hadrons at B- Factories 3 

varies from 0.31 pb at m c =1.2 GeV to 24 fb at m c =1.8 GeV. The cross section 
decreases as m c increases mainly because available phase space shrinks. If one can 
increase the cm. energy of the e + e~ , the m c dependence will decrease. 

In Fig. IHb) we show the differential cross section with respect to the invari- 
ant mass of cc. This is the prediction for da(e + e~ — > cc + X)/dm cc in leading 
order in a s . Experimentally, this differential cross section can be compared with 
the J2h H' da(e + e~ — > HH' + X)/dm,HH'i where H and H 1 are charm hadrons, 
which do not include anticharm. 

With a[e + e~ — > cc+X] ss 0.1 pb and current integrated luminocity C « 300 fb -1 
we expect roughly 30 events will be detected by the Belle detector. If our leading- 
order prediction is comparable to the measured value, it is very probable that the 
QCD higher-order corrections to the J/ip + r/ c cross section is small. Then the large 
discrepancy in J/tp + r\ c cross section may be due to the violation of factorization 
or existence of new production mechanism. If the measured cross section for the 
four charm hadron inclusive production is much larger than our prediction like the 
case of J/tp + r/ c , it is very likely that perturbative QCD corrections to J/tp + r\ c 
cross section is large enough to explain the discrepancy, which leads to the failure 
of reliability in perturbative expansion. 
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